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Motivated by the desire to construct meson-meson operators of definite relative momentum in 
order to study resonances in lattice QCD, we present a set of single-meson interpolating fields at non- 
zero momentum that respect the reduced symmetry of a cubic lattice in a finite cubic volume. These 
operators follow from the subduction of operators of definite heficity into irreducible representations 
of the appropriate little groups. We show their effectiveness in explicit computations where we find 
that the spectrum of states interpolated by these operators is close to diagonal in heficity, admitting 
a description in terms of single-meson states of identified J^ . The variationally determined optimal 
superpositions of the operators for each state give rapid relaxation in Euclidean time to that state, 
ideal for the construction of meson-meson operators and for the evaluation of matrix elements at 
finite momentum. 
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I. INTRODUCTION 

Hadron spectroscopy is concerned with the study of 
resonances which appear as poles of hadron scattering 
amplitudes. These scattering amplitudes are ultimately 
determined by the underlying theory of quarks and glu- 
ons known as QCD. Our best current general-purpose 
technique for computation of quantities in QCD is lat- 
tice QCD, which considers the field theory numerically 
on a grid of Euclidean space-time points of finite extent. 
Scattering amplitudes are not directly accessible in the 
Euclidean theory, but computations of the discrete spec- 
trum of states in a finite volume offer a way to access 
scattering amplitudes through a formalism presented by 
Liischer|T]. Within this framework, narrow meson res- 
onances appear through admixtures of localised single- 
hadron states, which are interpolated well from the vac- 
uum by local qq-\\ke operators, and meson-meson states 
having definite relative momentum determined by the 
finite-volume boundary conditions. 

In order to accurately extract a complete finite- volume 
meson spectrum, it is necessary to include in a basis of 
interpolating fields not just local operators resembling 
single-hadron states, but also some which resemble pairs 
of mesons projected into definite relative momentum. 
Our desire then is to form operators that can interpolate 
the meson-meson pair Ma Mb which will be of the form 
Oma{p)^Mb{p')7 where Oma{p) efficiently interpolates 
meson Ma with momentum p (and similarly Omb{p') 
efficiently interpolates meson Mb)- An optimal opera- 
tor for Ma with momentum p follows from variational 
diagonalisation of a matrix of single-hadron correlators 
evaluated with a basis of suitable operators. The subse- 
quent reduction of the contribution of excited Af^ mesons 
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to the MaMb correlators will mean that the relevant en- 
ergy levels, and hence signals for scattering, can be ex- 
tracted at earlier Euclidean times where statistical noise 
is smaller. 

This suggests that we need to develop a basis of single- 
hadron interpolating fields having non-zero momentum. 
In lattice QCD this cannot be achieved by simply Lorentz 
boosting operators at rest, as one is required to take into 
account the reduced symmetry of a discrete grid of finite 
extent. The lattice discretisation breaks rotational sym- 
metry at small distances, while the finite volume breaks 
rotational symmetry at large distances. The latter is of 
fundamental importance for states at non-zero momen- 
tum because the allowed momenta are determined by the 
boundary conditions which implement the symmetry of 
the "box". The relationship between the discrete finite- 
volume energy spectrum and the infinite volume scatter- 
ing amplitudes is also determined by the symmetry of the 
boundary. In most lattice QCD calculations, including 
those presented in this paper, both the lattice discreti- 
sation and the box have the same cubic symmetry, but 
this need not necessarily be so. For mesons in fiight, the 
symmetry is further reduced to the little group of allowed 
cubic rotations that leave the meson momentum invari- 
ant and it is operators transforming irreducibly under 
this reduced symmetry group that we will derive. 

As well as their use in constructing meson-meson op- 
erators, in-flight meson interpolators play a significant 
role in computations of meson matrix elements of certain 
currents, such as those used to measure electromagnetic 
and weak form-factors. Including operators which inter- 
polate mesons at non-zero momentum increases the num- 
ber of kinematic points available, allowing the form fac- 
tor to be calculated at more Q^ values. These quantities 
can reveal information about the quark-gluon structure 
of hadrons [H |3] . 

We will build on recent success [4] in constructing a 
basis of meson interpolators at rest that transform irre- 



ducibly under the cubic symmetry of the lattice. It was 
found that on the dynamical lattices of Refs. [5l[6], the 
effect of the finite lattice spacing was relatively small; 
operators suhduced from definite continuum spin into ir- 
reducible representations of the octahedral group over- 
lapped with states in a manner compatible with an effec- 
tively restored rotational symmetry. The observed spec- 
trum of states was explained in terms of single-meson 
states of determined J^'~^ with J^ values subduced across 
cubic irreps, and, as expected for a calculation featuring 
only fermion-bilinear operators, there was no clear obser- 
vation of a spectrum of meson-meson states whose dis- 
tribution across irreps is determined by the symmetry of 
the finite-volume boundary conditions. 

In this article we follow a similar philosophy, construct- 
ing fermion-bilinear operators for mesons in flight which, 
having a clear interpretation in an infinite volume con- 
tinuum, allow the continuum quantum numbers of single- 
hadron states to be identified. We then demonstrate the 
effectiveness of these constructions by using them to ex- 
tract an extensive spectrum of isovector mesons in flight 
at various momenta, with the continuum spins and par- 
ities reliably identified. With the extracted spectrum 
viewed in the finest detail we will observe effects that 
may be compatible with the admixture of meson-meson 
states. We propose that by including meson-meson oper- 
ators constructed from products of the in-flight operators 
derived herein, in future calculations we will be able to re- 
solve this admixture precisely and extract meson-meson 
scattering information. While Ref. [1^ considered only 
the case of the total system at rest, it has been shown 
that more information can be gained from calculations 
that work in a moving frame [71 [5]. 

We begin in Section In] by considering an infinite vol- 
ume continuum, constructing helicity operators and dis- 
cussing operator-state overlaps. In Section |III| we move 
to a finite volume, showing how these operators sub- 
duce into lattice little group irreducible representations 
to give subduced helicity operators. After giving a brief 
description of the lattice data sets used and the computa- 
tional method in Section [TV] in Section|V]we describe the 
method for identifying the J^'-^ of extracted states and 
show some examples of its application. The extracted 
spin-identified spectra of mesons in flight are given in 
Section IVll and we conclude in Section IVIIIl 



II. MESONS IN FLIGHT IN AN INFINITE 
VOLUME 

A. Helicity operators for mesons in flight 

We begin by describing the construction of helicity op- 
erators which enable the energies of mesons in flight to 
be extracted in an infinite volume continuum, before con- 
sidering the overlap of these operators with states. In 



duced into the irreducible representations relevant for a 
cubic lattice in a finite cubic volume, and how, because 
they have a definite infinite volume interpretation, they 
enable the identification of the spin and parity of states. 
As a starting point, consider the fermion-bilinear op- 
erators constructed in Ref. [1] , referring to that reference 
for more details, 



O'^^'ip)^ J2 CGs(mi,m2,m3, 

mi ,?7i2,?n3,... 



'E' 



^ip-x 



_ J V, J V, 

X V'(:?,t)rmi Dm^ Z^mj ...'0(f,t) . (1) 

Here D = D — D is a, gauge-covariant derivative, T is 
any product of Dirac gamma matrices, tjj is a quark field 
(smeared using distillation in our implementation) and 
spin, flavour and colour indices have been suppressed for 
clarity. As described in [5], the vector- like gamma ma- 
trices and derivatives are expressed in a circular basis 
and then coupled together using standard SU{2) Cleb- 
sch Gordan coefficients (represented by "CGs" in Eq. fTl) 
to form the operator O'^'^^p). However, for consistency 
between the operator constructions below and conven- 
tional definitions for helicity states, here we use a slightly 
different circular basis compared to that reference, 
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and similarly for the vector-like gamma matrices. Here 
ei{0,m) is the polarisation vector, given in Appendix pi 
for a spin-one meson at rest with spin z-component m. 

With full continuum rotational symmetry and p = 0, 
these operators have definite spin (J), spin ^-component 
(M) and parity, P (determined by the choice of F and the 
number of derivatives), i.e. they only overlap with states 
having these quantum numbers. In addition, if tp and tp 
correspond to quarks of the same flavour, these operators 
have definite charge conjugation parity C (determined 
by the choice of F, the number of derivatives and how 
these derivatives are coupled together), generalising to 
G-parity as appropriate. In the continuum the overlap of 
these operators with states is given by 



{0\O-''''{p^ 6)\p^ 0; J', Af) = Z^-'^Sjj,S_ 



M.AI' 



(3) 



Section III we will show how these operators can be sub- 



where the conserved parity and charge conjugation parity 
quantum numbers have been suppressed. 

At non-zero momentum, the spin z-component M is 
not a good quantum number unless the momentum is di- 
rected along the z-axis; in general, Eq. |3] is not applica- 
ble to py^ 0. It is more convenient to consider operators 
with definite helicity, the projection of the spin compo- 
nent along the direction of p. In Appendix |A] we show 
that, in analogy to the transformation from basis states 
with definite spin z-component to helicity states, helicity 



operators can be constructed by 
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(p) = E^ma(^)^''''(p1 
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where O"^' is a helicity operator with helicity A, 'D\j^{R) 
is a Wigner-I? matrix, and R is the (active) transforma- 
tion that rotates (0, 0, |p|) to p. Note that in general 
there are an infinite number of R which rotate (0,0, \p\) 
to plj There is some subtlety in the appropriate choice of 
R when the symmetry is reduced (e.g. on a finite volume 
lattice with a finite lattice spacing) and we discuss this 
below in Section ITlTl 



An equivalent way to construct these operators would 
be to choose an initial circular basis for the vector- 
like gamma matrices and derivatives defined in terms 
of the component of spin along p instead of the spin 
z-component. These vectors could then be coupled to- 
gether using the standard SU{2) Clebsch Gordan co- 
efficients as described above. This emphasises that by 
constructing helicity operators we've just effected a basis 
change. As will be illustrated below, the resulting basis is 
much more convenient for studying mesons with non-zero 
momenta because it respects the symmetries of a system 
with a meson in flight and helicity is a good quantum 
number (in an infinite volume continuum). 



B. Helicity operator overlaps 

Here we discuss the overlap of helicity states 
onto the helicity operators constructed above, 
<^0|O'^''^''^(p)|p; J'^ , A'), still considering an infinite 
volume continuum. 

The overlaps for operators with p = Q have a simple 
form, Eq. [sj J' = J, P' = P and A' = A. These con- 
straints arise from the rotational symmetry in 3 spatial 
dimensions of a system containing a particle at rest, and 
no further constraints are gained by imposing Lorentz 
symmetrjFl This can be seen explicitly by comparing the 
Lorentz covariant parameterisations given in Appendix 
A of Ref. [TT] and the 3-rotation covariant parameterisa- 
tions given in AppendixJD] herein. However, at finite mo- 
mentum, 3-rotation symmetry, a subgroup of full Lorentz 
symmetry, provides less stringent constraints, something 
that can again be seen explicitly by comparing those two 



^ If {<j>, 6, i/i) represents a rotation around the z-sxis by ■ip followed 
by a rotation around the y-axis by 9 and finally a rotation around 
the 2-axis by </>, one convention is ij = (</>, 8, 0) |9] and another, 
the Jacob-Wick convention 1101 , is i? = {(f), 9, —</>). Including an 
additional arbitrary initial rotation ip will still give a rotation 
from (0,0, IpI) to p. Different conventions will lead to different 
phases in the definition of states and operators. 

^ Although Lorentz symmetry constraints would relate operators 
containing the temporal and spatial pieces of Lorentz vectors, 
e.g. 7" and 7'. 



different parameterisations. Since in the operator con- 
structions we treat space and time asymmetrically, we 
will consider in the detail the restricted 3-rotational sym- 
metry constraints. 

Using only the constraints arising from 3-rotation sym- 
metry, the overlap of a state of definite J^ {P, P' are the 
parities at rest) and helicity onto a helicity operator is 
given by 



(o| 
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The helicity of the state. A', is constrained to be the same 
of that of the operator, A, but an operator constructed 
to have integer spin J at rest can in general overlap onto 
states of any integer spin J' when boosted to non-zero 
momentum. This is because at non-zero momentum the 
3-dimensional rotational group is broken to a subgroup, 
called the little group, made up of the rotations and re- 
fiections which leave the momentum, p, invariant. In an 
infinite volume continuum the little group is U(l) |12II13| . 
the group of rotations and reflections in two dimensions, 
regardless of the momentum direction. The irreducible 
representations of this little group are labelled by the 
magnitude of helicity (and also, for A = 0, a 'parity', 
?7, described below): lAI** = 0+, 0", i, 1, f , 2, |,3, . . . . 
Apart from the two one-dimensional irreps with A = 0, 
these irreps are all two-dimensional. Another conse- 
quence of this reduced symmetry is that an operator 
with |A| 7^ can overlap with states of both parities, 
and the overlap factors, Z, can depend on A. In Ap- 
pendix Id] we give the most general parameterisations of 
these operator-state overlaps and these features can be 
seen explicitly. Note that any flavour quantum number, 
such as charge conjugation parity or its generalisation, is 
still a good quantum number if it is so at rest. 

Helicity states at non-zero momentum are not eigen- 
states of parity because a parity transformation, P, re- 
verses the direction of the momentum, p — > —p. A reflec- 
tion in a plane containing the momentum direction, H 
(a parity transformation followed by a rotation to bring 
the momentum direction back to the original direction), 
preserves p. However, helicity states are also in general 
not eigenstates of H because under such a transformation 
A — > — A. If we consider a state with p — pi along the 
z-axis and a reflection in the yz plane {x — > — x), H 
then we have (Appendix [A 1| 

where J^ are the quantum numbers of the state at rest 
and fj = P(— 1)"^. It can be seen that for A = the 
helicity states are eigenstates of H with eigenvalue fj and 
therefore|3 

(0|O-^-^-^=°(p)|p;J'^',A') = 



yz^ 



z 



.1,J',P,P' 



'',A=0]^^ 



)),))'OA',0 



(6) 



^ and because A = there is no additional phase depending on the 
choice of reflection plane (Appendix I A 1\ 



For A 7^ the helicity states are not eigenstates of n q| 

If we had started with Lorcntz 4- vector gamma matri- 
ces and derivatives instead of 3-vectors, we could have 
constructed operators which only overlap onto states 
with one J^ at non-zero momentum. Performing a 
Lorentz boost on the state and operator on the left-hand 
side of Eq. |3] would leave the right-hand side unchanged, 
or putting it slightly differently, we can always boost 
back to the rest frame where Eq. [3] is valid. However, 
our use of only spatial derivatives precludes us from do- 
ing this. We do not include temporal derivatives in our 
operator constructions because in distillation the quark 
fields in the operators are only smeared in the spatial 
directions. This means that we can not form a lattice- 
discretised temporal derivative which is simply a rotation 
of a lattice-discretised spatial derivative. In addition, we 
have used anisotropic lattices and smeared differently in 
the temporal direction in the action, and these lead to 
similar problems. If we naively included temporal deriva- 
tives, these would be on a different footing to the spatial 
derivatives and so the resulting operators would not be 
Lorentz covariant and would overlap onto more than one 
J^ at non-zero momentum. 

Although our operators (and hence the operator-state 
overlaps) are not Lorentz covariant because they have 
been built from 3-vectors, in a theory with Lorcntz sym- 
metry there are still more stringent constraints on the 
overlaps than those given in Eq. [5] In Appendix [D] we 
show the restrictions on state-operator overlaps which 
arise from the constraints of a Lorentz symmetric theory 
applied to our helicity operators. The J^ of the states 
with which the operator has non-zero overlap depends 
on which gamma matrices and derivatives the operator 
is constructed from and how these have been coupled to- 
gether. An operator which at rest overlaps with only one 
J can in general overlap with states of many J at non- 
zero momentum, but the set of allowed J^ is reduced by 
Lorentz symmetry constraints compared to the infinite 
set allowed by 3-rotation symmetry, Eq. [5] 

To illustrate the above points, consider a simple ex- 
ample where the fermion bilinear operator consists of a 
spatial gamma matrix in the Cartesian basis and there 
are no derivatives. Projected onto zero momentum, this 
operator, Oi = i/'TiV'? has non-zero overlap with only the 
vector state, J^'~^ = 1 , 

(0|Oi(p = 0)|p = 0;l— (Af)) = Ze,(p-0,Af) , (7) 

where Si is the polarisation vector of the state and M is 
its spin z-component. Now consider a non-zero momen- 
tum along the positive z-axitrl pz = (0, 0,p). In a theory 



with only 3-rotation symmetry and not full Lorentz sym- 
metry, the state-operator overlaps are given in the last 
column of Table pQ| in Appendix |DJ i.e. the operator has 
non-zero overlap with states of all integer J . For states 
with J = 0, 1 the overlaps arer] 



(o|o,(p;)|p-p;;o+-) 

{0\Oi{pz)\p = Pz-l—{M)) 



ZnPi , 

Zle^ + Z[{e^Pj)p, 



However, in a theory with Lorentz symmetry most of 
these overlaps are constrained to be zero. For example, 
there are no Lorentz covariant structures that reduce to 
the Z[ [{e^pj)pi\ or Z^ [eijtP^e'^] terms above. In this 
case the operator has non-zero overlap with only two 
states (shown in the column labelled "LI" in the afore- 
mentioned table), 

{0\OM)\p = p'z;0+-)<xp^ , 

(0|O,(p;)|p-p;;l— (A/)) cxe,(p,M) . 

If instead we consider the Lorentz covariant operator 
Op = i^jfiip, the non-zero overlaps are 

(o|Op(p;)|p = p;;o+-)(xpp, 

(0|Op(p;)|p = p;;l— (M)) (xep(p,Af) . 

In this case we can take linear combinations of the four 
components to project onto just one of the states: tak- 
ing the inner product with p'* projects onto the 0^ 
state, whereas taking the inner product with e*^{p,M') 
projects onto the M' = M component of the 1 state 
only. 

In summary, although simplified somewhat by using 
helicity operators, the pattern of operator-state overlaps 
at non-zero momentum is more complicated than at zero 
momentum. In particular, because we construct opera- 
tors out of 3-vectors rather than Lorentz 4-vectors, we 
can not in general form operators that have non-zero 
overlap with states of only one J^ . However, as we shall 
show in Section Ivj the remaining constraints prove to be 
enough for us to be able to identify the states' spins and 
parities. In the following section we consider complica- 
tions arising from a finite cubic volume. 



III. MESONS IN FLIGHT AND HELICITY 
OPERATORS ON THE LATTICE 

In general, the symmetry of the lattice discretisation 
(which breaks rotational symmetry at small distances) 
need not be the same as the symmetry of the finite vol- 
ume, the boundary conditions (which break rotational 
symmetry at large distances). However, in this work we 



* although we can form eigenstates by taking hnear combinations, 
~ \p; J, A) it rjlp; J, —A), they are of Umited use here. 

^ chosen to be along the 2-axis so that the helicity state/operator 
is the same Jz state/operator 



note that the overlap onto the M = piece of the I"' 
a consequence of Eq. [6| 



vanishes, 
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helicity operator, from a helicity operator: 



TABLE I. Allowed lattice momenta on a cubic lattice in a 
finite cubic box, along with the corresponding little groups 
(the double covers relevant for integer and half-integer spin) 
from Ref . |12l 113] . We list only the single cover irreps rel- 
evant for integer spin. Lattice momenta are given in units 
of 2TT/{Lsas) where n,m,p £ Z* are non-zero integers with 
n ^ m j^ p. The A and B irreps have dimension one, E two 
and T three. Did is the dicyclic group of order 4n. 



only consider a cubic lattice in a finite cubic box with pe- 
riodic boundary conditions, and so both the lattice and 
the boundary have the same symmetry. In this case, 
the full rotational symmetry of the continuum is bro- 
ken to the (double cover) of the octahedral group, O^, 
equivalent to the (double cover) of the symmetry group 
of the cube. Spin is no longer a good quantum num- 
ber and states are not classified by J but instead by 
the irreducible representations, irreps (A), of O^. In 
Rcfs. m [T3], it was found that on the dynamical lat- 
tices of Refs. O |6] , the effect of the finite lattice spacing 
was relatively small. Operators with definite continuum 
spin (Eq. [I]) were subduced into octahedral group irreps; 
the observed spectrum of states and operator-state over- 
laps were compatible with an effective restoration of ro- 
tational symmetry, and the spectrum was explained in 
terms of single-meson states of determined J^'^ with J^ 
values subduced across these irreps. 

In going from zero momentum to non-zero momentum 
the symmetry is broken further. Whereas, in an infinite 
volume continuum the little group is the same regardless 
of the momentum direction, p, in a finite volume the par- 
ticular lattice little group depends on the star of p 12] . 
The allowed p are determined by the boundary conditions 
and the star of p is the set of all p related by allowed lat- 
tice rotations; as a shorthand, we refer to these different 
stars as different momentum types. Momenta are quan- 
tised by the periodic boundary conditions on the cubic 
box and we give all momenta in units of r^^ ; Os and 
Ls are respectively the spatial lattice spacing and spa- 
tial extent of the lattice in lattice units. The relevant 
little groups are given in Refs. [HI [13]. We summarise 
the allowed lattice momenta and the corresponding little 
groups and irreps in Table |l| 

In analogy to Ref. |31, we consider subduction coeffi- 
cients, 5JJ' , which specify how the helicity, A, subduces 
into a little group irrep A (row fi = 1 . . . dim [A]). Using 
these we can construct a little group operator, a subduced 



t^'^Hp) 



A=±|A| 



t> 



^J,P,X 



iP), 



(8) 



where tj — P(— 1)'^ with J and P the spin and parity of 
the operator 0'^'^'^{p — 0). The subduced helicity op- 
erators are different orthogonal combinations of the two 
signs of helicity, -|-|A| and — |A|. These subduction coeffi- 
cients can be calculated using the group theoretic projec- 
tion formula (for example, see Appendix A of Ref. ^). 

In Table |ll] we give these subduction coefficient^ for 
momenta of the form (n, 0,0), {n,n,0) and {n,n,n) 
[n S Z*) for |A| < 4. More details of our conventions 
are given in Appendix O For all these momenta, A'^ = 
0+ and 0~ subduce onto the Ai and A2 irreps respec- 
tively. However, the other A contained in those irreps 
(i.e. which A can mix with A = 0) depends on the 
momentum type: (n, 0,0) (Dic4) also has A — ±4,...; 
(n, n,0) (Dic2) has A = ±2, ±4,...; {n,n,n) (Dica) has 
A — ±3, .... Note that, although the pattern of subduc- 
tions does not depend on conventions, the relative phases 
in subduction coefficients can be convention dependent, 
determined by the helicity operator construction, the ro- 
tations used (discussed below) and the particular repre- 
sentation matrices chosen for the two-dimensional irreps. 

Instead of constructing little group operators via helic- 
ity operators which are then subduced into little group 
irreps, we could have projected directly from our contin- 
uum operators (in the J^ basis) into these little group 
irreps, but the subduction coefficients would then in gen- 
eral depend on the momentum direction and not just the 
momentum type. Alternatively, we could have worked 
only in terms of lattice irreps, subducing from octahe- 
dral group irreps to lattice little group irreps. However, 
the method we have described makes it clear how the 
operators can be physically interpreted in the limit of an 
infinite volume continuum and, as we show below, en- 
ables the continuum quantum numbers of states to be 
identified. 

As noted above, in a finite volume there is some sub- 
tlety in the choice of the R in Eq. H (out of an infi- 
nite number of possibilities) which rotates (0, 0, \p\) (not 
necessarily an allowed lattice momentum) to p (an al- 
lowed lattice momentum). In an infinite volume con- 
tinuum the choice is not important as long as one con- 
vention is chosen because the two directions perpendic- 
ular to p are interchangeable; they are related by a ro- 
tation. However, on a finite lattice the two transverse 
directions are not in general interchangeable; they are 
not related by an allowed lattice rotation. Consider, for 



^ Note that each subduction coefficient could be multiplied by an 
arbitrary phase and it would still give a subduction from the 
helicity to the little group irrep. 
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4 


^2(1) 


(5,,+ -775,,_)/V2 


Dies 


0+ 


^l(l) 


1 


(n, n, n) 


0~ 


^2(1) 


1 




1 


E2{\) 


(5,,+ ±»75,,_)/x/2 




2 


E2{1) 


(±(5,,+ -f?<5.,_)/x/2 




3 


^l(l) 


{Ss^+-vSs.-)/V2 




3 


^2(1) 


{Ss,+ +'n5s,-)/V2 




4 


E2{\) 


{5,,+ TvSs.^)/V2 



TABLE II. Subduction coefficients, 5^'j;, for |A| < 4 with 
s = sign(A); other notation is defined in the text. 



example, the two (shortest) lattice vectors perpendicular 
to (0,1,1), namely (1,0,0) and (0,1,-1) - these have 
different lengths. Therefore, if the i?'s are not chosen 
appropriately there is the possibility of an inconsistency 
between correlators from different momentum directions 
within the same momentum type. For example, if R used 
for (0, 1, 1) and R' used for (1, 1, 0) are not related by a 
lattice rotation, this could lead to effectively different 
definitions of the little group irreps and different bases 
for httle group rows in (0, 1, 1) compared to (1, 1, 0); e.g. 
Bi (B2) from (0, 1, 1) could correspond to B2 (-Bi) from 
(1, 1,0). This would prove problematic if averaging over 
different momentum directions to increase statistics and 
when constructing multi-meson operators which contain 
a sum over momentum directions with various weights. 
If R and R' are related by an allowed lattice rotation, 
these inconsistencies do not arise because there is no al- 
lowed lattice rotation that rotates from direction (1, 0, 0) 
to direction (0, 1, —1), and so the different irreps or irrep 
rows cannot be mixed up. 

We ensure consistency between different momentum 
directions by breaking down R into two stages: R — 
R\atRici- First, i?refi rotatc from (0,0, \p\) to Prof, where 
Prof is a reference direction for momenta of type p (i.e. for 



the star of p) . Note that i?rof is not in general an allowed 
lattice rotation, but this is permissible because all we are 
effecting by this rotation is a basis transformation. For 
-Riat we choose a lattice rotation which rotates from prof 
to p. There are in general a finite number of possible 
choices of lattice rotations for i?iat; the particular choice 
is not important as long as we make the choice consis- 
tently for each momentum direction (for example when 
later constructing multi- meson operators). More details 
of our implementation are given in Appendix [Cj 

To illustrate this, consider momentum directions with 
p^ = 2 and choose pj-of = (0, 1, 1)- If we consider p = 
(0, 1, 1) then R = R^ef is just a rotation from (0, 0,p) to 
p. On the other hand, if we consider p = (1, 1,0) then 
R = -Riat^rcf where i?iof is that just described and i?iat 
is a lattice rotation that rotates (0, 1, 1) to (1, 1, 0). 

Now that we have described the construction of sub- 
duced helicity operators, in the next section we briefly 
describe some of the lattice and computational details, 
before explaining our method for determining the contin- 
uum spin of extracted states and giving some examples 
in the following section. 



IV. COMPUTATIONAL DETAILS 

We use anisotropic dynamical Clover lattices with lat- 
tice parameters described in detail in Refs. [5l[6], having 
a spatial lattice spacing, as ^ 0.12 fm, and a tempo- 
ral lattice spacing approximately f = as /at — 3.5 tim es 

we 



III 



smaller, corresponding to a^ ^ 5.6 GeV. In Table 
give details of the data sets used; in all cases we increase 
statistics by averaging over equivalent momentum direc- 
tions within a given momentum type. In the main body 
of results shown here, we use the Nf = 3 data set with 
lattice size L^ x T = 16'^ x 128 in lattice units. This en- 
semble has three degenerate dynamical quark flavours, 
i.e. has SU(3) flavour symmetry, and corresponds to 
TO^ = niK = TTiri ~ 700 MeV. Some comparisons are per- 
formed using the lattices of different volumes and quark 
masses given in Table |III[ We consider isovector mesons 
and so only connected Wick contractions contribute. 

Correlators are constructed using the distillation 
method [15] which, in combination with these anisotropic 
lattices, has proved fruitful in a number of applications 
[H [T¥[ [TB1 - [TO| . For studying mesons in flight, an impor- 
tant aspect of distillation is that it allows a large basis 
of operators at the source and the sink, each projected 
onto a definite momentum. 

Our operator basis consists of all possible combina- 
tions of gamma matrices and zero, one or two lattice- 
discretised gauge-covariant derivatives, coupled together 
and subduced into lattice little group irreps, to form sub- 
duced helicity operators as described above. The nota- 
tion follows that of Ref. [4]: an operator (F x DjJ)'l 
contains a gamma matrix, F, with the naming scheme 
given in Table |IV[ and A^ derivatives coupled to spin Jjj , 
altogether coupled to spin J; A refers to the helicity com- 
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-0.0743 
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-0.0743 


-0.0743 
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(1,1,0) 


6(12) 
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-0.0743 


-0.0743 
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(1,1,1) 


8(8) 
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-0.0743 


-0.0743 
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16^ X 128 


64 


(2,0,0) 


6(6) 


200 


3 


-0.0743 


-0.0743 
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20^ X 128 
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(1,0,0) 
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-0.0840 


-0.0743 
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20^ X 128 
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(1,0,0) 


6(6) 
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6 


-0.0840 


-0.0743 
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20^ X 128 
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(1,1,0) 


6(12) 
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6 


-0.0840 


-0.0743 
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24^ X 128 
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(1,0,0) 


6(6) 
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4 


-0.0840 


-0.0743 


396 


24^ X 128 
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(1,1,0) 


6(12) 


553 


4 



TABLE III. The data sets used, rrii is the bare mass of the degenerate up and down quarks, and rUs is the bare strange quark 
mass. A^vccs is the number of distillation vectors used, A'^mom is the number of momentum directions averaged over out of a 
possible Ai'possibio mom, A'^cfgs is the number of configurations and Nt^^^, the number of time-sources. 



flO TT 712 bo p p2 fll 



6i 



r 1 75 7o75 70 7i 7i7o 757» (^ijk'yj^k 
TABLE IV. Gamma matrix naming scheme. 



ponent. 

The number of operators in each irrep is given in Table 
\V\ For the two-dimensional irreps we increase statistics 
by averaging the correlators over the two irrep rows. We 
consider isovector mesons which are eigenstates of charge 
conjugation (or G-parity) and so the irreps are labelled by 
charge conjugation parity, C, as well as little group irrep, 
A. Correlators were analysed using our implementation 
of the variational method[20 l [2T ] described in [4j. We use 
the fl baryon mass to set the scale and quote all energies 
as ratios, E/mci, where atmn = 0.353(3) on the lattices 
with m^ w 700 MeV ^. 

As a first test of the subduced helicity operator con- 
structions, for each possible momentum with |pj^ — 1, 2, 
3 and 4, all possible correlators were computed, including 
those off-diagonal in irrep and/or irrep row. This allowed 
the orthogonality of correlators between different irreps, 
and between different rows in the same irrep, to be veri- 
fied, along with the positivity of diagonal correlators and, 
within each irrep, the hermiticity of the correlator matrix 
and the consistency of different rows in the same irrep. 
In addition, the consistency of correlators from different 
momenta directions with the same momentum type was 
verified, something which would not necessarily be true 
if the rotation matrices, R, had not been chosen appro- 
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23 




BtW 
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At{l) 


21 




Et{-2) 


33 




^r(i) 


21 




^2^(1) 


15 




E^m 


35 



priately (see Section III | 



TABLE V. The number of operators in our basis for each 
little group irrep A'-' (C is the charge conjugation parity); 
also shown is the dimension of each irrep. 



V. SPIN DETERMINATION 

For mesons with non-zero momentum, the identifica- 
tion of the continuum spin of a state extracted in a lat- 
tice calculation is more complicated than for mesons at 
rest; even in an infinite volume continuum, an operator 
in general can have non-zero overlap with states of many 



J , as discussed in Section 



II B 



The reduced symmetry 



compared to the zero momentum case means that more 
states appear in each irrep leading to many degenera- 
cies in the spectrum; some are dynamical degeneracies of 
physical states in the true spectrum and others are caused 
by the reduced symmetry. Trying to identify the contin- 
uum quantum numbers by matching degenerate levels in 
different little group irreps (after taking the continuum 



limit) will be even more problematic than for mesons at 
rest. It would appear difficult or impossible to disentan- 
gle the states shown in Figs. [8| to [TT] without using some 
more information beyond solely their energies. 

For mesons stable against hadronic decay, a general 
idea as to where to expect different energy levels can be 
obtained from using the mass of the mesons at rest (ex- 
tracted on the same lattice) and the dispersion relation, 



{a,Em? = {cL,mf + 



1 /27r 



e \L 



bt 



(9) 



To determine a precise value for ^ = cis/cit on this lat- 
tice, we use our measured ir energies, ajE'^dpl), at |p|^ = 
1,2,3,4 and the tt mass from Ref. [TT, and fit the dis- 
persion relatiorFlto obtain a good fit with ^ = 3.441(8). 
However, because of the degeneracies (within statistical 
uncertainties) in the spectrum, some further informa- 
tion is needed to determine the spin of extracted states. 
Therefore, following a similar procedure to that devel- 
oped for mesons at rest in Refs. [Hm], we consider the 
state-operator overlaps. 

After subducing to lattice little group irreps, and ne- 
glecting the effects of a finite volume and finite lattice 
spacing, Eq. [5] becomes 



(o| 



.[J,P,|A|].^l ^ J, 



IP' 



\') = z 



.LJ ,P,P- 



'■''^sl 



A,p.^K^' 



(10) 



where fj = P(— l)'^ and J and P are respectively the 
spin and parity of the operator at rest. As discussed 
above, even in an infinite volume continuum, at non-zero 
momentum an operator can in general overlap onto states 
of many different spins and both parities. We use the 
constraints on operator-state overlaps which arise from 
Lorentz symmetry applied to our helicity operators to 
determine the continuum spin and parity of states. These 
constraints can be read off from columns "LI", "L2" and 
"L3" of Tables pOTll pQVl |XV] and pTyl] in Appendix |D] 

We observe that mixings due to the reduced symme- 
try are small; a subduced operator from a helicity opera- 
tor with helicity magnitude |A| only overlaps significantly 
onto states of helicity magnitude |A| and not with other 
helicities in that irrep. In Fig. IT] we show the normalised 
correlation matrix on time-slice 5 for the A^ irrep with 
IpI^ = 2 (including zero, one and two-derivative opera- 
tors) ordered such that those subduced from A = come 
before those from |A| = 2. The correlation matrix is ob- 
served to be approximately block diagonal, correlators 
mixing A = and |A| = 2 are small. 

This lack of mixing between different | A | is also seen in 
the Z values. For example, consider the A"-^ = Aj irrep. 
In Fig. [2] we show the Z's for the lightest five states in 
this irrep with |p|^ — 1, along with the J^^ assignment 



Fitting to a more general form of the dispersion relation with p^ 
and p^ terms gives a consistent result, g = 3.441(13), and does 
not improve the goodness of fit. 



A=0 




Al =2 



0.2 0.4 



0.6 0.8 1.0 



FIG. 1. Normalised correlation matrix (Cijj \JCiiCjj) on 
time-slice 5 in the A^ irrep with |pl^ = 2 (Dic2). The first 
20 operators are subduced from A = (top and left), the re- 
maining 11 operators are subduced from |A| = 2 (bottom and 
right). 



for each state; for simplicity of presentation, we use a 
smaller operator basis, only including the seven zero and 
one-derivative operators. The operator naming scheme 
is described in Section [TV] Also shown on the right of the 
figure is the spectrum of energy levels extracted in this 
irrep (the spin-identified spectrum is shown in the left 
hand panel of Fig. [s]). 

The same operators are shown for each state but the 
colour coding varies according to the J^'^ hypothesis: 
red bars correspond to operators which overlap onto 
states with that J^'^ at rest, blue bars correspond to op- 
erators which overlap with such states only at non-zero 
momentum and hatched grey bars correspond to oper- 
ators which do not overlap with such states (using the 
constraints arising from Lorentz symmetry). As a con- 
crete example, consider the first two operators shown, 
(7r)3[=Q and {t^i)'\^- These contain, respectively, 75 and 
7o75 in spin space, both with no derivatives, and at zero 
momentum both overlap only with J^'~^ = '" states. 
Both contain no vector indices (i.e. no 7' or D*) and so 
Table [Xni] is relevant but, as explained in Appendix [D| 
because these operators overlap with J^ = 0~ at rest, we 
must flip all the state parities in that table. The operator 
{ttYxZq is of the form given in the "LI" column and from 
there we see that, using Lorentz covariance constraints, 
it only overlaps onto J^ = 0^. In the figures it is there- 
fore coloured red for 0~ states and grey for other J^ . On 
the other hand, {t^2)'\Zo ^^ *-'^ ^^^ form given in the "L2" 
column and we see that at non-zero momentum it can 
overlap with both 0~ and 1+. It is therefore coloured 
red for 0~ states, blue for 1+ states and grey for other 

For all the states shown, it can be seen that the red 



bars are largest, although these overlaps can be small if 
there is a mismatch in structure between the operator 
and the state. The blue bars can be significant, these 
should be proportional to some power of \'p\/M and so 
are suppressed compared to the red bars at low momen- 
tum, and the grey bars are small. From these figures it 
can be seen that in each case the J^'-^ can be determined 
unambiguously. For |p|^ = 1 (Dic4) the continuum helic- 
ities which subduce into the A2 irrep are 0, 4, . . . (Table 
In]) . The restricted operator basis used here does not con- 
tain helicities greater than two and so we can not test 
hclicity mixing here. In any case, the lightest spin-four 
state is expected to be significantly higher in energy than 
the heaviest states we extract in this analysis. 

In Fig. [3]we show the Z^s for the overlap of the lightest 
six states in the A J irrep with |p|^ = 2 onto the ten zero 
and one-derivative operators in that irrep, along with 
the spin assignment hypothesis for each state. Again, for 
each state, the continuum J^*^ can be unambiguously 
determined. For |p|^ — 2 (Dic2) the continuum helici- 
ties which subduce into the A2 irrep are 0, 2, 4, . . . (Table 
|ll[ ) and the ^2 operator basis contains operators coming 
from A = (the first seven operators) and |A| = 2 (the 
last three operators). It can be seen that for each state 
only A = or I A| = 2 overlaps are statistically signif- 
icantly non-zero, further supporting the observation at 
the level of the correlation matrix. Fig. [l] Furthermore, 
these overlaps are generally smaller than the overlaps 
onto operators with the correct continuum helicity but 
which are constrained to be zero by Lorentz symmetry. 

A second stage of the spin identification procedure is 
to compare the Z's for the overlap of a given state onto 
a given continuum operator when this operator is sub- 
duced into two or more different little group irreps. From 



Eq. 10 neglecting any finite volume or discretisation ef- 
fects, these overlap factors should be equal, something 
that is empirically found to be true for mesons at rest 
extracted on this lattice in Refs. [11[T3]. For mesons in 
flight the opportunities for such comparisons are more 
limited. The only directly analogous situation is when 
the same |A| subduces into two one-dimensional irreps 
(as opposed to subducing into one two-dimensional ir- 
rep); for example, A = ±2 with |p|2 = 1 subduces into 
Bi and B2, and A = ±1 with \p\^ — 2 subduces into Bi 
and B2. These comparisons can only be used to confirm 
the continuum helicity rather than the continuum spin 
(although determining the continuum helicity does give 
a minimum value of the continuum spin) and so such 
comparisons are less useful for mesons in flight than for 
mesons at rest. 

As an example, in Fig. |4]we show the Z's for the low- 
est two J^*-^ = 2=*=+ states with |pj^ = 1 overlapping 
onto the three one-derivative operators subduced into the 
A'-' = B^ and i?^ irreps. The good agreement between 
Z values with the same operator subduced into differ- 
ent irreps is apparent. In Fig. [5] we show the Z's for 
the lowest two J^'^ = 1*~ states with |p|^ = 2 over- 



lapping onto the eight zero and one-derivative operators 
subduced into the A*-^ = B^ and i?^ irreps. Again, there 
is good agreement between the Z's in different irreps. 

The consistency of these Z values for the same |A| sub- 
duced into different little group irreps relied on three- 
rotation symmetry, assuming that the effects of a flnite 
volume and flnite lattice spacing are small for these quan- 
tities. Using only constraints arising from three-rotation 
symmetry, the Z's for different |A| components overlap- 
ping onto the same continuum operator are not related, 
but some of these can be related using Lorentz symme- 
try constrai nt^ These constraints can be read off from 
Tables pOnl pOVl pCV] and pCVll in Appendix [D) Because 
the constraints depend on which gamma matrices and 
how many derivatives the operator was constructed out 
of, and how these were coupled together, comparisons 
are much more involved compared to those for mesons 
at rest. Therefore, in general we do not consider these 
comparisons as part of our spin-identiflcation procedure. 

As an example, consider operators containing a 7* 
and zero derivatives, along with operators containing one 

derivative (-D*) and the identity in spin-space. Reading 
off from Table pQV ("LI" column), we see that the over- 
lap onto a J^ = 1~ state is iZi for A — ±1 and i^Zi for 
A = 0. Note that if we only imposed 3-rotation covari- 
ance, the A = overlap has two independent constants, 
Zi and Z2, and so cannot be related to the A = ±1 over- 
laps. In Fig. [6] we show the overlap factor for both these 
operators for the lowest 1 state with |p|^ = 1. There 
is good agreement between the different irreps for the 7* 

operator (1); for the D^ operator (2) the Z values are 
similar but disagree signiflcantly outside their statistical 
uncertainties. Sources for this discrepancy could be vi- 
olations of hypercubic symmetry arising from the action 
or the result of mixing with multi-meson states which 
have a different distribution across irreps compared to 
single-meson states. To check whether this discrepancy 
could be a flnite volume effect, in Fig. [6] we also plot the 
Z-values scaled by •y/16'^/20'^ for the corresponding state 
extracted on a lattice with the same quark masses but a 
larger volume, 20^ in lattice units. Although the Z values 
for operator (2) get slightly closer together on the larger 
volume, there is no signiflcant change in the pattern. 

In summary, we have shown how operator overlaps can 
be used to identify the continuum J^'~'' of mesons with 
non-zero momentum extracted in lattice calculations. In 
Appendix [E| we show some tests of the robustness of our 
extracted spectra and in the following section we apply 
this method to identify continuum spin in our full ex- 
tracted spectra. 



^ if our operators were built out of Lorentz vectors we could relate 
all the helicities in every case, but our use of operators built from 
3- vectors limits where we can relate different |A| 
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FIG. 2. Overlaps, jZj, of operators onto states labelled by E/mu in the A — A2 irrep with \p\'^ — 1 (Dic4); the extracted 
spectrum is shown on the right. The Z's are normalised so that the largest value for that operator across all states is equal to 
1. The unshaded area at the head of each bar represents the one sigma uncertainty on either side of the mean. For each state 
the bars represent, from left to right, the operators (7r)^zg, {■K2)i=o, (ai)A=o> (P x D''}Li)i=o, (Pa x D^jU)i^h, {h x Dy^lj^=g 
and (61 X I^j1i)a=o- The colour coding is explained in the text. 
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FIG. 3. As Fig. M but for the A — A^ irrep with |p]^ — 2 (Dic2). For each state the bars represent, from left to right, 
the operators (7r)^=g, {n^YxZl, {a^YxZl, (p x D^l^)izl, (p2 x djliYx^l (^1 x D^l{)izl, {h x D^jU)i=l {p x D^jU)(^%, 
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FIG. 4. \Z\ values for the A = ±2 components of J^'^ = 
2="=+ states with |p|^ = 1 (Dic4) subduced into the A^ = B+ 
and B2 irreps. The inset shows the extracted lowest-lying 
energy levels. The operators are: (1) (p x Dj1j)|^T£2i (2) 
(P2 X D^jl^)]'^^^ and (3) (fei x D^}1^)(^^2; the Z's from the 
Bf irrep (shifted to the left) and the B^ (shifted to the right) 
are plotted. 
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FIG. 5. \Z\ values for the A 

states with |p|^ = 2 (Dic2) subduced into the A*^ — B^ and 
B2 irreps. The inset shows the extracted lowest-lying en- 
ergy levels. The operators are: (1) {p){x^Li, (2) {p2){x\Li, (3) 
{bi)(,^L^, (4) {n X Df^,){,^l„ (5) (^2 X D^jU){,^U, (6) (ao x 
D^}U)wU, (7) (ai X D^jlXt, and (8) {a, x D^jU)(,%; 
the Z's from the B^ irrep (shifted to the left) and the B^ 
(shifted to the right) are plotted. 



VI. SPECTRA OF MESONS IN FLIGHT 



Using the data sets given in Table III with m. 



700 

MeV and volume 16'^ x 128, and the operator basis de- 
scribed above (including zero, one and two derivatives, 
and enumerated in Table lv| , for each irrep and momen- 
tum type we computed a matrix of correlators. These 
correlators were analysed using the implementation of 
the variational method [IHl HI] described in Ref. [4 and 
we identified continuum spins using the methodology de- 



scribed above. 

Before we show the full extracted spectra for each |p|^, 
we illustrate some features by considering one case in 
more detail, the lowest lying states with positive charge 
conjugation parity for |p|^ — 1 (Dic4), shown in Fig. u\ 
The colour coding indicates the identified J^: black/grey 
for J = 0, red for J = 1, green for J = 2 and blue foiFj 



there are no J = 3 states in Fig. [t] but these appear in later 
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FIG. 6. \Z\ values for the J^*^ — 1 ground state with 
IpI^ = 1 (Dic4) subduced into the A'^ = A^ {X = 0) and 
E2 (A = ±1) irreps. The operators are: (1) (p)''^^ ~ 7', 
(2) (ao X dJ/Ii)-'=^ ~ D'; the Z's from the A^ irrep and Z's 
scaled by E/M from the E2 irrep are plotted. For comparison, 
also shown are the Z values, all scaled by (■yT&^~/20^), for a 
larger lattice with volume 20'' in lattice units. 



J = 3; darker shades with a solid outline correspond 
to positive parity and lighter shades with a dashed out- 
line correspond to negative parity. Orange boxes with 
a dotted outline indicate that the state's J^ could not 
be unambiguously determined in this analysis. The lines 
in the right hand column show expected energies (with- 
out statistical uncertainty) for the lower-lying states us- 
ing masses from Ref. JTl] and the dispersion relation, 
Eq. [9) with ^ = 3.441 (Section |v]). The triangles (cir- 
cles) show the distribution of the components of these 
expected states with definite J^=+ (J^^~) across the 
irreps, using Table [11} 

Spin assignment was initially carried out indepen- 
dently for each irrep and, as can be seen, the different 
components of the same continuum state can be matched 
up between different irreps with good agreement between 
the energies. The pattern of states across irreps (the 
boxes in the figure) is the same as the pattern expected 
for states of definite continuum J^ subduced across ir- 
reps (the triangles and circles in the figure) . We note that 
it would be extremely challenging to identify the different 
states with E /m^ ^ 1.34 if we did not use information 
from the Z values. 

In the A^ irrep the state with E/nin ^ 0.90 (dark 
grey, solid outline) is identified as having J^'-' — 0++, 
has no partner in the other irreps and lies where the 
flo is expected. The next state in this irrep ~ 1.01 (dark 
green, solid outline) is identified as the zero helicity piece 
of 2+"^p^ It can be matched up with states identified as 
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2++ in the E} (A = ±1) and the B+ and B+ (A = ±2) 
irreps, and lies where the 02 is expected. The state at 
^ 1.29 (orange, dotted outline) is not well determined 
in this analysis. There are two states ~ 1.34, one (dark 
grey, solid outline) is identified as O"*""*" and appears to be 
the Oq. The other (light red, dashed outline) is identified 
as the zero helicity piece of the exotic 1 ^ and can be 
matched up with the 1 '" in the E2 irrep (A = ±1). This 
exotic state is extracted with the same precision as the 
non-exotic states; its considerable overlap onto operators 
proportional to the commutator of two covariant deriva- 
tives, requiring a non-trivial gluonic field configuration, 
suggests a hybrid interpretation [H [32] . 

Moving to the Aj irrep, there are three states (light 
grey, dashed outline) which have been identified as 
JPC ^ Q-+ at E/mn - 0.53, 1.08, 1.34; these do not 
have partners in other irreps and lie where the pion and 
excited pious are expected. The state at ^ 0.97 (dark red, 
solid outline) is identified as the zero helicity piece of l"*"^, 
matches up with the 1+"*" in the E2 irrep (A = ±1), and 
lies where the Oi is expected. Finally, the state ~ 1.25 
(light green, dashed outline) is identified as the zero helic- 
ity piece of 2 '", can be matched up with the 2 ^ states 
in the E^ (A = ±1) and the B^ and B^ (A = ±2) irreps, 
and lies where the tt2 is expected. 

The full spectra with |pj^ = 1,2,3 and 4 are shown in 
Figs. [8j [9j [To] and [TT] respectively, and the patterns can 
be understood in a similar way. The expected distribu- 
tion of the components of states with definite continuum 
jP across the little group irreps can be determined from 
Table ]n] and the extracted states are observed to match 
this pattern. 

We are able to extract states with exotic quantum 
numbers, identified as hybrid mesons. In addition, the 
third excited 1 state {p'"), E/nin ^ 1-4 at |p|^ = 1, 
has considerable overlap onto operators proportional to 
the commutator of two covariant derivatives, requiring a 
non-trivial gluonic field configuration. This suggests that 
we identify this state as a non-exotic hybrid mesorl^and 
corresponds to the non-exotic hybrid identified at zero 
momentum [31 122). 

The reduced symmetry of the Dic2 and Dic3 little 
groups compared to Dic4 leads to even more degeneracies 
in the spectrum. For example, in the left hand panel of 
Fig. p\ there are two degenerate states in A^ irrep with 
E/mn ~ 1-07 identified as J^^ = 2++; one is the A = 
piece and one is a A = ±2 piece, identified from their Z 
values. The other A — ±2 piece is in the A^ irrep and the 
A = ±1 pieces are in the Bf and S^ irreps. A similar 
situation arises in the A2 irrep with two states ~ 1.29 
identified as 2 *", one corresponding to A = and the 
other to A = ±2. 



figures 
^^ Recall that the A = piece of a state with fj = P(— I)'' = +1 or 
— 1 subduces into the Ai or A2 irrep respectively, but the |A| > 
pieces of a state subduce into the same irreps regardless of parity. 



^^ or a mixture of hybrid and conventional mesons with a significant 
hybrid component 
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We note that, around E/ttiq ~ 1.4 there start to be 
states that we can not reliably extract or identify in some 
irreps. For example, in A2 and A^ irreps with |p|^ = 1, 
we are missing the A = pieces of the 3+"*" and 3"' states 
respectively (Fig. [s]). This could be because, although 
we have an extensive operator basis, there is only one 
operator which overlaps with J = 3 at rest in each of 
those two irreps. To reliably extract a spin-three state 
and increase the number of other states we can extract, 
we would need to extend our operator basis to include 
three-derivative operators. 



In Fig. 12 we show the extracted energies as a function 
of IpI^ for some of the lower- lying states. Also shown are 
the expected energies given by the dispersion relation, 
Eq. |9j using ^ = 3.441 and the rest mass extracted in 
Ref. j4J . It can be seen that the extracted energies are in 
reasonably good agreement with the dispersion relation, 
highlighting the observation that the extracted energies 
are in good agreement with the expected states in Figs.[8J 
|9J [To] and IT and indicating that with the dynamical 
tuning of parameters in the anisotropic lattice action, 
symmetry between space and time has been restored to 
good accuracy. 

In summary, for all the |p|^ considered, we are able to 
identify the J^*^ of an extensive spectrum of extracted 
states across all irreps. The observed pattern of states 
across irreps is the same as the pattern expected for the 
components of states with definite continuum J dis- 
tributed across irreps (Table In]). The energies of the dif- 
ferent helicity components can be matched up between 



irreps and generally agree well. However, the compo- 
nents in different irreps do not have to be degenerate in 
a finite volume and a possible interpretation for this is in 
terms of mixing with multi-particle states, something we 
discuss in the following section. 



VII. MULTI-MESON STATES AND HELICITIES 

In an infinite volume continuum with rotational sym- 
metry in three dimensions, a state's two helicity compo- 
nents with the same |A|, +|A| and — |A|, are constrained 
to have equal energy; they correspond to two rows of 
the same two-dimensional irrep. However, components 
with different magnitudes of helicity, |A|, of a state are 
not constrained to have equal energy; they correspond to 
different irreps. It is the additional constraints imposed 
by Lorentz symmetry, discussed in Section |IIB[ which 
require that all the different helicity components have 
the same energy; it is in this sense that they correspond 
to one state, a state at rest which has been boosted to 
non-zero momentum. 

In a finite volume, the reduced symmetry means that 
there are even fewer constraints and, depending on the 
little group corresponding to the particular momentum 
type, the two helicity components -f|A| and — |A| may 
be split in energy. For example, with |p|^ — 1 (Dic4), 
the two A = ±1 components both subduce into the two- 
dimensional E2 irrep and so are constrained to have equal 
energies. However, with |p|^ = 2 (Dic2), the two A = ±1 
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components are split between the B\ and B2 irreps and 
so are not constrained to be degenerate. 

In the spectra we have extracted, in general we see very 
little splitting of energies between the different irreps. As 



an example, in the left hand panel of Fig. 13 we show the 
lowest energies extracted in the A~^ and E'j" irreps with 
IpI^ = 1, and in the A^, B^ and i?^ irreps with |p|^ = 2, 
along with the expected energies for the p [J^'-^ — 1 ) 
in flight from the dispersion relation. There is no sta- 
tistically significant splitting between different |A| or be- 
tween the B^ and i?^ energy levels, and these energies 
lie where expected from the dispersion relation. We note 
that the lowest energy tttt states with overall momenta 
p= (0,0,1) andp= (0,1,1) are 7r(0,0,0)7r(0,0, 1) and 
7r(0, 0, 0)7r(0, 1, 1) respectively. In the absence of interac- 
tions, these would have energies E/mi-i w 0.95 and 1.04 
respectively, significantly heavier than the corresponding 
ground state energies we extract in the relevant irreps. 
This suggests that any mixing with tttt states will be 
small. In fact, this is one reason for testing the helic- 
ity operator formalism on lattices with a relatively heavy 
pion. 

For comparison, in the right hand panel of Fig. [13] 
we show the corresponding energies extracted on 20^^^ 
and 24^^ x 128 lattices with lighter up and down quarks 
{Nf = 2 + 1), corresponding to mjr ~ 400 MeV and hav- 
ing atmn = 0.2951(22) [15] and ^ = 3.433(5) (from a fit 
to the TT dispersion relations across all volumes rj Here 
we observe more significant splittings between the ener- 
gies in A^ compared to E2 or i3j~,_B^, and a smaller 
splitting between Bi and B2] the pattern is the same on 
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both volumes. These splittings could arise from an im- 
perfect restoration of Lorentz symmetry; however, com- 
parison between the two lattices suggests an explanation 
in terms of mixing with multi-meson states in a finite 
volume. 

On the lattice with smaller TOtt and volume 20'^ [24^], 
the lowest energy tttt state with overall momentum p — 
(0, 0, 1) is 7r(0, 0, 0)7r(0, 0, 1) and in the absence of inter- 
actions this has energy E/rnn sa 0.63 [0.59]; the low- 
est TTTT state with overall momentum p = (0,1,1) is 
7r(0, 0, 0)7r(0, 1, 1) at « 0.73 [0.67] in the absence of inter- 
actions. The relative closeness of these tttt energy levels 
to the energies we extract suggests that the tttt admix- 
ture in will be larger here compared to the lattices with 
heavier pions. 

For p = (0,0,1), in the absence of interactions the 
7r(0, 0,0)7r(0,0, 1) has lowest energy {E/rrifi w 0.63 or 
0.59 on the 20'^ or 24*^ respectively) and it subduces 
only into the Ai irrep. The next lowest energy state 
is 7r(0, -1, 0)7r(0, 1, 1) (« 0.89 or 0.79) and this subduces 
into both the Ai and E2 irrepj^ Therefore, the pattern 
of TTTT energy levels in the Ai irrep is different from the 
pattern in the E2 irrep and so the effect of multi-meson 
mixing in the p will be different in the two irreps leading 
to a splitting of energies. 

A similar interpretation applies to p = (0,1,1): the 
state 7r(0, 0, 0)7r(0, 1, 1) has lowest energy [E/ma ~ 0.73 
or 0.67) in the absence of interactions and it subduces 
only into the Ai irrep. The next lowest energy state is 
7r(0, l,0)7r(0,0, 1) (« 0.78 or 0.70), subducing into Ai 
and Bi, and the next is 7r(0, l,0)7r(l, 1, 1) (w 0.98 or 
0.86), subducing into Ai and -62- Therefore, we can again 
qualitatively explain the splitting of the energy levels in 
terms of mixing with multi-meson states. 

In summary, multi-particle states have a different pat- 
tern of energy levels across irreps compared to a single- 
particle state. The pattern of degeneracies seen in Figs.jS) 
|9j [To] and [TT] do not show any clear evidence for multi- 
particle states. However, there are suggestions that as we 
reduce the pion mass mixing with such states will become 
more important. 

In the full spectrum, we should observe additional 
states, the orthogonal mixture of single and multi-meson 
basis states. However, we do not expect to be able to 
properly resolve multi-particle states with only fermion 
bilinear operators. As discussed in Ref. [1] , if the opera- 
tor basis is not sufficient, the variational method can not 
separate the states and a conservative approach is to sug- 
gest that our mass values are only accurate up to the or- 
der of the hadronic width of the state. By supplementing 
our operator basis with two-meson operators (built from 
the helicity operators described herein), we can extract 
the full spectrum of single and two-meson states. We will 



^^ For the analyses on the lattices with m-n- ^ 400 MeV, the opera- 
tor basis consisted of only zero and one-derivative operators. 



^* and also the B^ irrep; this will be expanded upon in future work, 
but for the present discussion it suffices to observe that the pat- 
tern of subductions into Ai is different from E2 ■ 
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then be able to extract phase shifts and, if appropriate, 
resonance parameters using the method of LiischerjT], 
generalised to an overall non-zero momentum "7", utilised 
explicitly in ^; this will be the subject of future work. 



VIII. CONCLUSIONS 

We have shown how to construct fermion-bilinear op- 
erators for mesons in flight which are subduced from 
operators of definite helicity and respect the reduced 
symmetry of a cubic lattice in a finite cubic volume. 
We have demonstrated the effectiveness of these oper- 
ator constructions by using them to extract an extensive 
spectrum of isovectors mesons with various momenta. 
We find that the spectrum can be described in terms of 
single-meson states of determined J^ and have reliably 
identified the quantum numbers of the extracted states, 
including states with exotic quantum numbers, states of 
high spin and excited states. 

In order to study resonances and scattering in lattice 
QCD, it is necessary to extract a complete finite-volume 
meson spectrum, and to do this we must supplement 
our bases of operators resembling single-hadron states 
with operators resembling pairs of mesons projected onto 
definite relative momentum. Such operators, interpolat- 
ing the meson-meson pair Ma Mb , will be of the form 
Oma{p)^Mb{p'): where Oma{p) efficiently interpolates 
meson M^ with momentum p. Using large bases of op- 
erators of the form described herein, we can variationally 
diagonalise the matrix of correlators and hence obtain the 
optimal operator for M^ with momentum p (and simi- 
larly for Mb)- The subsequent reduction of the contribu- 
tion of excited M^ mesons to the MaMb correlators will 
mean that the relevant energy levels can be extracted at 
earlier Euclidean times where statistical noise is smaller, 
leading to more precise determinations of signals for scat- 
tering. 

With the extracted spectrum at light pion masses 



viewed in the finest detail, we observe effects that may be 
compatible with the admixture of meson-meson states. 
In future calculations, meson-meson operators will be 
constructed from products of the in-flight operators de- 
rived herein, allowing us to resolve this admixture pre- 
cisely and so to extract information on meson-meson 
scattering and any possible resonance [U [3 |H| . 

As well as their application to constructing multi- 
meson operators, these in-flight interpolators play a sig- 
nificant role in calculations of matrix elements at finite 
momentum, such as those used to measure form fac- 
tors. The inclusion of meson operators with non-zero mo- 
mentum in such analyses gives further kinematic points, 
enabling the form factor to be calculated at more Q^ 
points and so the Q^ dependence to be mapped out more 
fully [H |3]. Again, having an optimal operator to inter- 
polate a meson will reduce the contribution of excited 
mesons, meaning that quantities can be extracted at ear- 
lier Euclidean times where statistical noise is smaller. 
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Appendix A: Helicity states and helicity operators 

Here we describe the construction of helicity operators^ 
using helicity states as an analogy and following the he- 
licity state constructions of Jacob and Wick [TD]. Parts 
of Chung [9 and Chapter 3 of Martin and Spearman [25] 
are also relevant. 

The operator R^^^^ = g-i04g-»eJ„g-»VJ. rotates 
around the z-axis by -0, then around the j/-axis by 9 and 
finally around the z-axis by (with a fixed coordinate 
system). As a shorthand we will denote the (active) ro- 
tation which rotates from the z-axis to the direction of p 
[which is {p,9,(/)) in spherical coordinates] by Rq. There 
are many possible choices for Rq and these lead to dif- 
ferent phases in the definition of helicity states. Jacob 
and Wick choose i?0,6i,<^ and Chung chooses R,f,.efl, but 
the choice is not important for the present discussion as 
long as one convention is used consistently. 

Denoting a boost along the z:-axis with magnitude p 
by Lzip), a helicity state is defined by 



\p;J,X)=RoL,{p)\J,X) , 



(Al) 



where | J, A) is a state at rest with spin J and J^- 
component A. From the definition, it is simple to show 
that under some arbitrary rotation, _R, these helicity 
states transform aa^ 

R\p; J, a) = \Rp; J, A) , 

and the helicity is invariant. 

In contrast, a canonical state (using the terminology 
of Chung), where the spin component is measured along 
the z-axis, is defined by 



\p; J, m) = RqLz{p)Rq ^ | J, m) 
= L{p)\J,m) , 



(A2) 



where L{p) = RqLz{p)Rq is a Lorentz boost along direc- 
tion p with magnitude p. These states transform under 
rotations as 

R\p; J,m)^Yl ^nJrniR)\m J, m') , 



where 'DI^/^^R) is a Wigner-2? matrix, i.e. these states 
transform under rotations in the same way as | J, m) . 

Using the definitions in Eqs. |Al| and |A2[ it is straight- 
forward to show that helicity and canonical states are 
related by 



|p;J,A)=^P^'l(i?o)|p;J, 



(A3) 



^^ Note that there may be some phase introduced here if the rota- 
tion has some component around the p direction. 
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We have shown previouslyfl^ that, in the cont inuum, 
when p = the O"^'*^ constructed in Section II A overlap 
with only one spin: 

(0|O^'"(p = 0)|p = 0; J',™') = Z^'^6j,j.6m,m' , (A4) 

and so we can write 

|0; J,to) = [O^'™(0)]^|0) . (A5) 

In particular, the operator, [C''''"'(0)]''', transforms in a 
similar way to the state. We can therefore carry over the 
discussion of canonical and helicity states and define a 
boosted canonical operator^ 

[0-'^™(p)]t ^ L(p)[O'^-'"(0)]tL-i(p) . (A6) 

A helicity operator is then given by 

P'-\p)]^ = E ^ml(^) [0''"(p1]^ , (A7) 

ni 

and transforms in the same way as a helicity state. Tak- 
ing the complex conjugate we obtain 

n'-'\p) ^J2'^l2*iR) ^'■"'ip) : (A8) 

rn 

which is Eq. [4] in Section II A 



1. Parity and reflections 

Helicity states at non-zero momentum are not eigen- 
states of parity because a parity transformation, P, re- 
verses the direction of the momentum, p — > —p. A reflec- 
tion in a plane containing the momentum direction, 11 
(a parity transformation followed by a rotation to bring 
the momentum direction back to the original direction), 
preserves p. However, helicity states are also in general 
not eigenstates of H because under such a transformation 

As an example, consider a helicity state with spin J, 
rest parity P and momentum, pi, directed along the z- 
axis, \p'z; J^, A); the state at rest is denoted by | J^, A). 
If Hyz is a reflection in the yz plane {x — >■ —x) and P is 
a parity transformation (x, y, z ^ —x, —y, —z) then 



H 



yz 



— i-rrj^p __ —iTTjy—i-Rj^p 



In analogy to the discussion in Ref. [TU], 

e-'^-^^e— -^^ |J^,A) =5]l?ifl(^,7r,0)|j^,A') , 

= (-l)'^|j^-A) . 

Since P commutes with R and L^ip) commutes with Tlyz, 
we have 



If instead we consider a reflection in the xz plane 
{y -^ — 2/), ^xz, then the discussion proceeds in a similar 
way [TU] . 



H., 



g-i7rj„p 



We have 



and this leads to 



Xz \pz; J'', A) = Pi-iy-^ {p-; JP, -A) 



tiyz IpI; J^: a) = P{-iy IpI; J^, -a) 



(A9) 



(AlO) 



This has the same form as Eq. |A9| but with a different 
phase. If we consider an arbitrary p and an arbitrary 
reflection in a plane containing p, the result will be of 
the same form as Eqs. |A9| and |A10| but in general there 
will be some A-dependent phase. 

It can be seen that for A = the helicity states are 
eigenstates of 11 with eigenvalue fj = P{—1)''. Further- 
more, because A = there is no additional phase depend- 
ing on the choice of reflection plane. The fj 'parity' is a 
good quantum number of the zero-hclicity components at 
flnitc momentum and determines which little group irrep 
A = subduces into (Table Ih]). In addition, because fj 
is related to the transformation of helicity states with 
A > under reflections, fj also appears in the subduction 
coefficients for A > 0. 



Appendix B: Helicity polarisation vectors and 
tensors 

Here we construct a representation for the polarisa- 
tion vectors of a spin-one helicity state and describe the 
generalisation to polarisation tensors. 

Taking the defining equation for e'(p, A) to be 

{0\O'\p;J=l,X)=Ze'{p,X) , (Bl) 

if follows that 

{0\R-^RO'R-^R\pz; J = 1, a) = Ze\pz,X) , (B2) 

and so, from the transformation properties of vectors and 
helicity states under rotations, 

{R-^y^{0\O^p;J=l,X)=Ze\pz,X) , (B3) 

and therefore 

{0\O'\p; J = 1, A) = ZR'^e^pz, A) . (B4) 

This implies that 

e\p,X)^R'^e'{Pz,X), (B5) 

and equivalently, since _R*^ is real, 

e"ip,X)^R''e^*iPz,X). (B6) 
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An explicit representation for the polarisation vectors 
of a spin-one helicity state with an arbitrary momentum 
can be constructed using Eqs. |B5| and |B6| and the stan- 
dard representation for p = pi, 



£'i{Pz,m^ 0) 
ei{pz,m = ±1) 



{0,0, E/M) , 
T^(1,±.0) 



(B7) 



where M and E are, respectively, the mass and energy 
of the state. 

Polarisation tensors for states of higher spin can be 
constructed in the usual way from direct products of the 
spin-one polarisation vectors with the helicities coupled 
together using standard SU{2) Clebsch-Gordan coeffi- 
cients. 



Appendix C: Lattice rotation and little group 
conventions 



In Table |VI| we give the specific rotations that we 
use for Rrci (Section III) which rotate from (0,0, \p\) to 
Prcf- We use the same convention as in Appendix \A\ 
namely that a rotation R^,B,ip — e^^'^'^^e^'^^'^^e^'^^'-'' ro- 
tates around the z-axis by ■(/', then around the y-axis by 
9 and finally around the z-axis by (j) (with a fixed coor- 
dinate system). 



In Tables |VII[ |VIII| and |IX| we give a choice of rep- 
resentation matrices for, respectively, the little groups 
Dic4, Dic2 and Dics. For Dic2 and Dica, the rotations 
and reflections refer to a coordinate system which has 
been transformed using i?rcf , so that p defines the new z- 
axis. These representations are used with the group the- 
oretic projection formula (for example, see Appendix A 
of Ref. [4 ) to calculate the subduction coefficients given 
in Table in 



Appendix D: Tables of general parameterisations 

Here we give general parameterisations for the infinite 
volume continuum overlaps of states with fermion bilin- 
ear operators using constraints of Lorentz or 3-rotation 
symmetry. Tables IX XI and XII give parameterisations 
for operators with zero, one and two vector indices re- 
spectively. The columns labelled "Lorentz" give Lorentz 
covariant parameterisations for the overlaps, i.e. us- 
ing Lorentz covariance constraints applied to our oper- 
ators built from 3-vectors. These are derived from the 
Lorentz covariant parameterisations given in Appendix 
A of Ref. [TT] by setting the Lorentz indices to specific 
values. The columns labelled "3-Rotation" give the gen- 
eral 3-rotation covariant parameterisations. 

In Tables XIII||XVT we give the overlap of states onto 
the operators O"'''^ defined in Eq. |4J These overlaps follow 
from Tables |X] |XI| and |XII| once the vector indices are 
transformed to a circular basis and then coupled together 



to give a definite spin as described in Section II A We 
only show overlaps where the helicity of the operator is 
equal to that of the state; the other overlaps are zero. 

The tables give overlaps for an operator of one parity, 
determined from the parity of the state with which it 
overlaps at zero momentum. To use the tables for opera- 
tors of the opposite parity, all the state parities must be 
fiipped; we give examples of this below. The inclusion of 
a 75 matrix in the operator does not change the allowed 
overlaps, except that the parity of the states may flip as 
discussed below. 

Note that the parameterisations are written in terms 
of helicity polarisation tensors, discussed in Appendix 
[B} not polarisation tensors defined in terms of the z- 
component of spin. Repeated indices are summed over. 
Although these are all fermion bilinear operators, for 
brevity we omit the ip and ip in the tables. 

As an example, consider the operator -ipj^ip which at 
zero momentum overlaps only onto the J^^ = 0^ state. 
The Lorentz covariant parameterisations are |llj 



(0|V'7^V|0+~(p)) ^ Zpf" 



and 



{0\^j^i^\l—{p,X))^Ze'^{p,X). 

Setting /i = these give 

(0|Vi7^=V|0+-(p)) =Zp°, 

and 

(0|V;7^=V|1""(P,A)) = Ze'>{p,\) = Zs\p,X)p'/E, 

where E — p'^ is the energy of the state and in the last 
line we have used p^e^(p. A) — 0. These two parameteri- 
sations are given in the "L2" column of Table pC| where we 
have absorbed factors of E into the arbitrary constants. 
As another example, consider the operator ipj^j^ip. 
This is the same as the operator considered previously 
apart from the additional 7^, and at zero momentum it 
overlaps only onto the J^'^ = '" state. In order to use 
Table |X] we must flip the parities of all the states and so 
the parameterisations become 



(0|V7^7^V|0~+(p)) - Zp", 



and 



{0\^j'-f^i;\l++{p,X)) = Ze'{p,X) = Ze\p,X)p'lE. 



Appendix E: Robustness of the extracted spectrum 

To check the robustness of our results to changing the 
operator basis, as well as using the full bases contain- 
ing up to two-derivative operators, we also performed 
some ana lyse s using only zero and one-derivative opera- 
tors. Fig. 14 shows the extracted energy levels in the A^ 
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tt] 1.4■ 



1.2■ 



1.0■ 



0.8 



0.6- 



At bf = 1 




0, 1 deriv 0,1.2 deriv 



0,1 deriv 0,1,2 deri 



FIG. 14. Comparison of the lower-lying states extracted us- 
ing zero and one-derivative operators, and with that basis 
supplemented with two-derivative operators, for the A2 irrep 
with IpI^ = 1 (Dic4) (left panel) and \p\^ = 2 (Dic2) (right 
panel). The box height shows the one sigma statistical uncer- 
tainty above and below the central value; the colour coding, 
indicating the J^, is described in the text. 



irrep for |p|^ = 1 and |p|^ — 2. All use 200 configurations, 
3 time-sources and the correlators are averaged over 6 
momentum directions, except for the one-derivative basis 
with \p\^ — 2 which is averaged over 12 momentum direc- 
tions. The one-derivative basis for this irrep contains 7 
(10) operators with |p|^ = 1 (|p|^ = 2), whereas the two- 
derivative basis contains 20 (31) operators. The colour 
coding indicates the identified J^: black/grey for J = 0, 
red for J = 1 and green for J — 2; darker shades with 
a solid outline correspond to positive parity and lighter 
shades with a dashed outline correspond to negative par- 
ity. Orange boxes with a dotted outline indicate that the 
state's J^ could not be unambiguously identified. 

Including two-derivative operators significantly in- 
creases the size of the operator basis, reduces the sta- 
tistical uncertainty on the higher extracted levels, and 
increases the number of states and the highest spin which 
can be reliably extracted (the latter point is not visible in 
these figures). However, for the lower states there is no 
statistically significant change in the extracted energies 
and no significant increase in precision of the energies 
extracted, verifying the robustness of the energies we ex- 
tract under a reduction in size of the operator basis. 
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Little Group 


Pref 





e 


V. 


DiC4 

Dic2 

Dies 


(0,0, n) 
(0,n,n) 

[n, n, n) 




^/2 
7r/4 




7r/4 
cos-i(l/V3) 




-^/2 




TABLE VL Rotations, -Rrcf, used, as described in the text. 



Irrep 


I 


R{7V) 


R(3n/2) 


Ri^m 


n 


i?(7r)n 


i?(7r/2)n 


R{3n/2)n 


Ai 


1 


1 


1 


1 


1 




1 


1 


A2 


1 


1 


1 


1 


-1 


-1 


-1 


-1 


-Ea 


(-1 0-1 

Vo i-* 


(-0^-^) 


(-°.T) 


iVo) 


il-\) 


Vol'' 


Cl'o) 


(-°.o) 


Bi 


1 


1 


-1 


-1 


1 




-1 


-1 


B2 


1 


1 


-1 


-1 


-1 


-1 


1 


1 



TABLE VII. Choice of representation matrices for the Dic4 little group. / denotes the identify transformation, R{(j>) denotes a 
rotation around the z-axis by and 11 denotes a reflection in the yz plane {x — > —x). 



Irrep 


I 


R{it) 


n 


i?(7r)n 


Ai 


1 


1 


1 


1 


A2 


1 


1 


-1 


-1 


Bi 


1 


-1 


1 


-1 


B2 


1 


-1 


-1 


1 



TABLE VIII. As Table |VII| but for the Dic2 little group. With these representations for B\ and B2 we are able to reproduce 
the multiplication table given in Ref. |13j . but we note that there is the possibility of another convention where Bi and B2 are 
interchanged. 



Irrep 


I 


R(2n/3) 


i?(-27r/3) 


R(n)n 


R{TT/3)n 


i?(57r/3)n 


A2 
E2 


1 
1 

ill) 


1 


1 

(-4^ 4 ) 


1 
-1 


1 
-1 

/ 1 ^■y3\ 
2 2 1 

\ 2 2 / 


1 



TABLE IX. As Table [yn] but for the Dies little group. 



State 


Lorentz 




3- Rot at ion 


J^ 


7 (LI) 


7" (L2) 




I or 7° 


0+ 


^1 


Zi 




^1 


1" 


- 


ZlEiPt 




ZlEiPr 


2+ 


- 


- 




ZxEii + Z2EijPiPj 


3- 


- 


- 




ZlEttjPj + Z2£tjkPiP]Pk 


4+ 


- 


- 


Zl £iijj 


+ Z2£zzjkPjPk + ZiEijklPiPjPkPl 



TABLE X. General Lorentz covariant and 3-rotation covariant parameterisations for the infinite volume continuum overlaps of 
operators with zero vector indices onto states up to spin-four. Zi are arbitrary constants, Pi is the 3-momentuni, for brevity 
we omit the arguments of the polarisation tensors e''''(p, A), and repeated indices are summed over. The parities given in 
the first column may all have an overall sign flip depending on which operator is considered: the parity of the state at rest 
overlapping with the operator must be determined and then the relative parities of the other states can be read off the table. 
More explanation is given in the text. 
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State 


O' (LI) 


Lorentz 

7"D" (L2) 


7°y (L3) 




3-Rotation 

0% -/°D' or 7%' 




0+ 

0- 


Z^p' 


^iP* 


- 


Zip' 


1+ 
1" 


Zie' 


Z±(eV±eV) 


Zie^'^PjEk 
Z-(e°p' - e'p°) 


Zie'^'^p.Ek 
Zie-' + Z2{e,pi)p' 


2+ 
2" 


- 


^le'^P, 


- 


Zie'^p, + Z2p'e^' + Zsp'e^'^p.p^ 
Zie'^'^PjEkiPi 


3+ 
3" 


- 


- 


- 


Zie''^''pj£kii + Z2e^^''pjekimPiPm 

Zip^e^'^Pk + Z2p'£^'"pjpm + ^3£*'pfePi + Zi£'^^ 


4+ 
4" 


- 


- 


- 


Zip'e""" 


+ Z2p'e"'''pkPi + Zap^e'^'^prnPiPm + ^4e""'=Pfc + Zse'' 

Zie'^''pjekttmPm + Z2("^''PjeUmnPlPmPn 


kl 

PjPkPl 



TABLE XL As Table |x| but for operators with one vector index. O' stands for 7' or D*. e'-''^ is the antisymmetric tensor. 
Where Zi± appears in the decomposition it means that both + and — terms are allowed. 



State 


O'' (LI) 


Lorentz 

-y^D'D' (L2) 


7°7'D^ (L3) 


3- Rot at ion 

0*^ -y^D'D' or j^fD' 


0+ 


ZiS'' + Z2p'p' 


Zi5'' + Z2p'p' 


Zi5'' 


ZiS'' + Z2PV 


0- 


- 


Zie'^'^Pk 


Zie'^'^Pk 


Zie^'^Pk 


1+ 


Z^e'^^ipkEo-poek) 


Zie''''p^Pkei + Z2e''''ek 
+Z±(e^'=ypfce,±6*^'=Pfee'p' 


Zi, Z2, Z- 

terms 
from "L2" 


Ziiie^'^pkEip' ± e'^'pkSip') + Z2€'"'ek + Z^e'^^pkEipi 


1" 


Zi±(eV±eV) 


Z3S'^e''pt: + Zip'e' 
+Z±{p°p°p^e'±pVs''pk) 


Z3,Z^ 

terms 
from "L2" 


Zi±{e'p' ± e^p') + Z2p'p'ekPk + Z^S'^e^Pk 


2+ 


Zie'^ 


Zse'^P^Pk 
+Z±(e^Vpfc±e''pV) 


Z- 

term 

from "L2" 


Zie'' + Z2pVe'"' + Zs^ie'^p^ ± e'''p')pk + Z^S'^s"" 
+Zr,S''e'"pm + Z.p'p^e'^'pm + Zre"''e^'""pkPmein 


2" 


- 


^ee'^Pfce^' 

+Z±(-6^'='p,e«pV 

±e'^'=e'='p,p«p" 

Te'^'pfce'>iPm) 


Zg,Z^ 

terms 

from "L2" 


Zi±{e''''pke''±e^'='pke") 

+^2±(e"''pfeP^e'>™ ± e^'='pfcpV'>™) 

+Z36"^"'=e'^'p, + Z4e'^''pke" + Z.^e'^^pke'^^pipr,, 


3+ 


- 


- 


- 


Zi±{e''='e''"' ± e'^'e^'^PhPrn + Z2±{e*'p^ ± e^'^'pOp^e'™ 

+^3±(e""V ± e^'='pOe'™>fcP™P„ + z^eiH"'' 
+Z,,e'^\'"^pW^ + Zee^'^Pke'^^^piPrnPn + Z^e'^^p^e'^^pi 


3- 


- 


Zie'^^'pk 


- 


Zie'^^Pk + Z2±(e"=V ± e^'V) + ^3±(e''V ± e'^'p^PhPi 

+ZipVe'='"^pmP^ + Z5p'p^e'"'pk + Z^S'^e^'^puPiPm 

+ZrS'^e'''W + ^86""e^™"e'"'-pfcP,nP,. 



TABLE XIL As Table |X| but for operators with two vector indices and states up to spin-three. O*-* stands for 'y'D-' or D^D-' 
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Operator 


state J^(A = 0) 


3- Rotation 


J (LI) 


7° (L2) 


(0,0,0 


0+ 


^1 


Zi 


Zi 


qO.O 


1" 


^1^ 


- 


Zi 


qO.O 


2+ 


/i^[Z.+i?^Z.] 


- 


- 


qO.O 


3" 


y^^[Z.+iJ^Z.] 


- 


- 


qO.o 


4+ 


2\/i ^ [Zi + E^Z2 + E^Zs] 


- 


- 



TABLE XIIL Overlaps of states (up to spin-four) witli the operators 0''°^'^ defined in Eq. l4l wiiere tiiese operators iiave been 
constructed from / or 7" and using tiie parameterisations in Table |x[ As in Table bcl "LI" and "L2" correspond to the Lorentz 
covariance constraints for operators constructed from, respectively, / and 7". The Zi in the "LI" and "L2" columns give the 
terms from the "3-Rotation" column which are allowed under these Lorentz constraints. Only the overlaps onto the zero helicity 
component of the states are shown, the overlaps with non-zero helicity are zero, p^ = \p\^ is the square of the 3- momentum, 
and M and E = p" are, respectively, the mass and energy of the state. 



Operator 


State J^ 








J- Rotation 






0^ (LI) 


7°I>" (L2) 


7"y (L3) 


0^'^ 


0+ 




0, ipZi, 




Zi 


Zi 


- 


0^'^ 


0^ 


- 


- 


- 


- 


0^'^ 


1+ 




pZi, 0, -pZi 




- 


Zi 


Zi 


0^'^ 


1" 




iZx, i§Zi+ii^Z2, iZ-i 




Zi 


Zl,Z2 


Zi = -E^Z2 


0^'^ 


2+ 


't&Z„ ^-^[E^Z,+p^Z,+p'E^Z,), ^z,] 


- 


Zi 


- 


0^'^ 


2- 


[1.0,-1] (^)^. 


- 


- 


- 


0^'^ 


3+ 


[1.0-1] (7&)(^i+^^^^) 


- 


- 


- 


O^'^ 


3- 


(J^)[^(£^^3 + ^4), 
ff (p2Zl-Ki?V^2+£^^3+^4), 


- 


- 


- 



TABLE XIV. Overlaps of states (up to spin-three) with the operators ([}-^°P'-^ defined in Eq. Ul where these operators have 
been constructed from operators with one vector index and using the parameterisations in Table pCl] "LI", "L2" and "L3" 
correspond to the Lorentz covariance constraints for the operators in Table [XI] the Zi in these columns refer to the terms in 
the "3-Rotation" column which are allowed under these Lorentz constraints. Only the overlaps where the helicity of the state 
is the same as that of the operator are shown; the others are zero. The basis is A = —1, 0, -f 1. 
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Operator 


state J^ 




3- Rotation 




O'' (LI) 


-y^D'D' (L2) 


7"yi3^ (L3) 


O0,o 


0+ 


^(3Zi+p2Z2) 


Zl,Z2 


Z\, Z2 


Zx 


0°'" 


1" 


^(2Zi++p2z2 + 3Z3) 


Zi+ 


Z\+, Z2, Z-i 


Z-A, 
1Z\-{- ^ —E Z2 


0)0,0 


2+ 


+3Z4 + 3^2 Z5 + E^p^Ze, - M^Zt] 


^1 


Z\ , ^3+ 


2Z3+E^ = -Zi 


O0,o 


3- 


y^f^ [^1 + 2Z2+ + 2E^Z3+ + E^p'Zi 
+p''Zs + 3E^Z(i + 3Z7 - M^Zsl 


- 


^1 


- 


O^'^ 


0+ 


- 


- 


- 


- 


o'-^ 


0" 


[0,iV2pZi,0] 


- 


^1 


Zi 


O^'^ 


1+ 


iV2 


'-{p'Zr- - ^2), f (^2 +p2^3),_(/2'i_ _ Z2) 




Z2 = -S'^3 


Zi_, ^2, -^3 


Z\-, Z2, Zz 


O^'^ 


1" 


V2pZi-[l,0,-l] 


Zx- 


Zx- 


Zi- 


O^'^ 


2+ 


t^^3- [1,0,-1] 


- 


Z,- 


Z,- 


O^'^ 


2" 


^[-E{Zr-+p^Z2--Z:,), 

^(-M^Zi- +p^Z4 + E^Z3 + E^p^Zs), 

-E{Z^^+p'Z2--Z3)\ 


- 


Z3 = -£^^5 


Za = -E'^Zt, 


O^'^ 


3+ 


f (-M^Zi- + Z4 + p'^Zj + S'Zs + E^p^Ze), 
- s/l{p^Z2- - Zi + E'p^Zz- + E^Z^- - £^^5)] 


- 


- 


- 


O^'^ 


3" 


2v^i^(^-+ ^^^3-) [1,0,-1] 


- 


- 


- 



TABLE XV. Overlaps of states (up to spin-three) with the operators 0''°^''^ defined in Eq. Ul where these operators have 
been constructed from operators with two vector indices and using the parameterisations in Table [Xll] "LI", "L2" and "L3" 
correspond to the Lorentz covariance constraints for the operators in Table [XIl] the Zi in these columns refer to the terms in 
the "3-Rotation" column which are allowed under these Lorentz constraints. Only the overlaps where the helicity of the state 
is the same as that of the operator are shown; the others are zero. The basis is A = — Jop, . . . , +Jop. Continued in TablelXVl] 
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Operator 


state J^ 




3-Rotation 




O'' (LI) 


j°D'D' (L2) 


-y^-y'D' (L3) 


O^'^ 


0+ 


[0,0, -sJ^p'^Zi, 0,0] 


Z2 


Z2 


- 


O^'^ 


0" 


- 


- 


- 


- 


O^'^ 


1+ 


iy^p'Zi+lO,!, 0,-1,0] 


- 


Zi+ 


Zi+ 


O^'^ 


1" 


-V2p 


0,Zi+, ^ {2Zi++p^Z2),Zi+,0 




Zi+ 


Zi+ , Z2 


2Zi-\- ^ —E Z2 


O^'^ 


2+ 


[-Zi+p^Z7,-§{Zi+p^Zs+), 

-3h{iM^+'2E^)Zi+2p^Z2 + 

4p^E^Z3+ + 2E^p''Z(, + M^p^^t), 

-§j {Zi+p^Z3+),~Zi+p^Z7] 


Zi 


Z\ , Z-i+ 


2Z3+E^ = -Zi 


O^'^ 


2- 


ip[2Zi+,§{Zi++p''Z2+), 

0,-§{Zi++p^Z2+),-2Zi+'j 


- 


Zi+ 


Zi+ 


O^'^ 


3+ 


%[eZ,+ , ^,j,{E^Z,+ +p^Z,+ + eVZ,+ }, 
0,-\/lJl{E^Zj++p^Z2+ + E^p^Z3+},-EZi+] 


- 


- 


- 


O"'^ 


3- 


-vki [e{Zi-P^Zs), ^^^{E'Z, +p^Z2+ + E^p'Z^+j , 
+2£;2p^Z4 + 2p^Z5 + M^'p^Zs}, 

^y/j{E^Zi+p^Z2+ + E''p'Z3+},EiZl -p^Zs)] 


- 


Zi 


- 



TABLE XVL Continuation of Table |X\^ 



